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Abstract. We study numerically the effect of varying the equation of state (EOS) on 
the formation of stellar clusters in turbulent molecular clouds. Our results show that the 
EOS plays an important role in the fragmentation of the clouds, and the determination 
of the initial mass function (IMF) of the protostellar cores. Fragmentation decreases 
with increasing adiabatic index 7 in the range 0.2 < 7 < 1.4, although the total amount 
of mass accreted appears to remain roughly constant through that range, resulting in 
more massive cores at higher 7. Fragmentation and collapse ceases entirely for 7 > 1.4 
as expected from analytic arguments. Primordial gas may have effective 7 > 1, in 
which case these results may help explain why models of the formation of the first 
stars produce isolated, massive objects. 

1 INTRODUCTION 

Fragmentation of molecular clouds is key to the formation of star clusters, yet 
the conditions for fragmentation are poorly understood. There have been a few 
analytical approaches (Jeans 1902; Penston 1966; Low & Lynden-Bell 1976), 
and numerical investigations of the effects of various physical processes on the 
collapse, such as geometry and rotation of the clouds (see Bonnell & Bastien 
1993 for a review), and magnetic fields (Galli et al. 2001). Recently, the effects 
of turbulence have been studied extensively in a series of 3D simulations (Klessen, 
Burkert, & Bate 1998, Paper I; Klessen, Heitsch, & Mac Low 2000, Paper II; 
Heitsch, Mac Low, & Klessen 2001, Paper III; and Klessen 2001, Paper IV). 
However, these results are based on isothermal models with 7=1. The effect 
of the equation of state (EOS) , which is essential in understanding the physical 
structure and stability of the turbulent clouds, remains uncertain. 

The balance of heating and cooling in a molecular cloud can be approximately 
described by the polytropic EOS, P = Kp" 1 , where if is a constant, and P, p and 
7 are thermal pressure, gas density and polytropic index, respectively. The den- 
sity structures generated by supersonic turbulence depends on the EOS (Scalo et 
al. 1998; Passot & Vazquez-Semadeni 1998; Passot & Vazquez- Semadeni 1999). 
A detailed analysis by Spaans & Silk (2000) suggests that 0.2 < 7 < 1.4 in the 
interstellar medium. 
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Fig. 1. (left) Number of sink particles representing collapsed star-forming cores as a 
function of time, and (right) fraction of total mass in the box accreted onto cores as a 
function of time, for different 7. 



In this paper, we perform a 3D numerical survey over this range of 7, to 
investigate the effects of the EOS on the fragmentation of turbulent clouds and 
the subsequent formation of a protostellar cluster. 



2 SIMULATIONS 

We carry out smoothed particle hydrodynamics (SPH) simulations to determine 
the effects of EOS on fragmentation. The code is implemented with periodic 
boundary conditions (Klessen 1997) and "sink" particles (Bate, Bonncll, & Price 
1995) to replace high-density cores. The driven turbulence is set up following the 
k =1-2 model in Mac Low (1999). For detailed descriptions of the code, see Paper 
II. The treatment of the polytropic cloud in the code is the same as in Bonnell 
& Bastien (1991). We vary 7 in steps of 0.1 in otherwise identical simulations. 
Self-gravity is included in the calculation at t = 2.0 after turbulence is fully 
established (see Paper II). Our models use 200,000 particles. The 13 simulations 
were performed simultaneously using a serial version of the code on processors 
of the Parallel Computing Facility of the AMNH. 



2.1 Fragmentation 

Figure 1 shows that the fragmentation properties and the accretion of the sinks 
strongly depend on the value of 7. The lower the 7, the earlier fragmentation 
occurs, and the more fragments form. It takes longer for gas with high 7 to 
fragment. The slope of the accretion curves are roughly the same, but the number 
of sinks is very different for different 7 at the same time, which suggests that 
individual protostellar cores is different for different 7. In low-7 gas, there are 
more cores, but the mass of each core is much smaller than in high-7 gas. 
The Jeans mass for a polytropic gas sphere Mj — (K-k/G) 272 p2 (7- 3) so 

dMj _ 37 - 4 (3 7 - 6 )/2 



dp 



OC (1) 
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Fig. 2. Top: 3-D distribution of the gas and protostars for different 7. Bottom: Mass 
spectra of gas clumps (thin lines) and of collapsed cores (hatched thick-lined histograms) 
for the corresponding cube above. The percentage shows the fraction of total mass 
accreted onto cores. The vertical line shows the SPH resolution limit. Shown also are 
two power-law spectra with v — —1.5 (dashed-line) and v = —2.33 (dotted line). 

which implies that 7 = 4/3 is a critical value during collapse. When 7 < 4/3, 
the Jeans mass decreases as density increases, promoting fragmentation, but 
when 7 > 4/3, the Jeans mass increases with density, stabilizing collapsing gas. 
Indeed, our model with 7 = 1.4 shows no collapse. 

2.2 The Mass Spectra 

Observations show that the mass distribution of gas clumps in various clouds 
follows a simple power law, dN/dM = M v , with v w —1.5 (Blitz 1993), while 
the stellar IMF follows v = —2.33 at the high-mass end (Salpeter 1955). 

Figure || shows that the spectra of both the clumps and the cores change 
with 7, but that the ones of cores change more dramatically. In low-7 models, 
the mass distribution of the cores at the high-mass end is roughly log-normal. As 
7 increases, fewer but more massive cores emerge. When 7 > 1.0, the distribution 
is dorminated by high mass cores only, and the spectrum tends to flatten out. 
It is no longer fitted by either a log-normal or a power-law. The clump mass 
spectra, on the other hand, do show power-law behavior at the high mass side, 
even for 7 > 1.0. 

Our results suggest that stars tend to form in clusters in a low-7 environment. 
Protostellar cores are low mass. The apparent lack of power-law behavior for the 
cores in the protostellar cluster might imply that simple accretion is unable to 
generate as many high-mass stars as predicted by the observations, which further 
suggests that other mechanisms such as collisions (Bonnell, Bate & Zinnecker 
1998) may be at work to produce the massive stars in a cluster. Higher resolution 
models will be necessary to confirm this, however. 
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On the other hand, our results also imply that massive stars can form in 
small groups or alone in gas with 7 > 1.0. Spaans & Silk (2000) suggest that a 
stiffer EOS (7 > 1.0) leads to a peaked IMF, biased toward massive stars, while 
an EOS with 7 < 1.0 results in a power-law IMF, in general agreement with our 
simulations. 

3 ISOLATED STARS 

The formation of isolated massive stars is of great interest, since in most cases, 
massive stars are found in clusters. But recently, Lamers et al. (2002) reported 
observations of isolated massive stars or very small groups of a few isolated mas- 
sive stars in the bulge of M51. Massey (2002) also reported many field massive 
stars in both the Large and Small Magellanic Clouds. From our simulations, we 
see that when 7 > 1, only very few or possibly only one fragment occurs. This 
cloud core is massive, and would result in the formation of a high-mass star. 

High resolution simulations by Abel, Bryan & Norman (2002) of the forma- 
tion of the first star suggest that initially only one massive metal-free star forms 
per pregalactic halo. In the early Universe, inefficient cooling due to the lack of 
metals may result in high 7. Our models would then suggest weak fragmentation, 
resulting in the formation of only one massive star per cloud. 

We thank M. Fall, H. Lamers, and H. Zinnecker for valuable discussions, and D. 
Janies and W. Wheeler for their work on the Parallel Computing Facility. YL thanks 
the AIP for its warm hospitality, and the Kade Foundation for support of her visit 
there. This research was partially funded by NASA grant NAG5-10103, and by NSF 
CAREER grant AST99-85392. 
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